Parallel Unsmoothed Aggregation Algebraic Multigrid Algorithms on GPUs

We design and implement a parallel algebraic multigrid method for isotropic graph Laplacian problems on multicore Graphical Processing Units (GPUs). The proposed AMG method is based on the aggregation framework. The setup phase of the algorithm uses a parallel maximal independent set algorithm in forming aggregates and the resulting coarse level hierarchy is then used in a K-cycle iteration solve phase with a $\ell^1$-Jacobi smoother. Numerical tests of a parallel implementation of the method for graphics processors are presented to demonstrate its effectiveness.Comment: 18 pages, 3 figure

A Jacobi-Davidson type method with a correction equation tailored for integral operators

The final publication is available at Springer via http://dx.doi.org/10.1007/s11075-012-9656-9We propose two iterative numerical methods for eigenvalue computations of large dimensional problems arising from finite approximations of integral operators, and describe their parallel implementation. A matrix representation of the problem on a space of moderate dimension, defined from an infinite dimensional one, is computed along with its eigenpairs. These are taken as initial approximations and iteratively refined, by means of a correction equation based on the reduced resolvent operator and performed on the moderate size space, to enhance their quality. Each refinement step requires the prolongation of the correction equation solution back to a higher dimensional space, defined from the infinite dimensional one. This approach is particularly adapted for the computation of eigenpair approximations of integral operators, where prolongation and restriction matrices can be easily built making a bridge between coarser and finer discretizations. We propose two methods that apply a Jacobi–Davidson like correction: Multipower Defect-Correction (MPDC), which uses a single-vector scheme, if the eigenvalues to refine are simple, and Rayleigh–Ritz Defect-Correction (RRDC), which is based on a projection onto an expanding subspace. Their main advantage lies in the fact that the correction equation is performed on a smaller space while for general solvers it is done on the higher dimensional one. We discuss implementation and parallelization details, using the PETSc and SLEPc packages. Also, numerical results on an astrophysics application, whose mathematical model involves a weakly singular integral operator, are presented.This work was partially supported by European Regional Development Fund through COMPETE, FCT-Fundacao para a Ciencia e a Tecnologia through CMUP-Centro de Matematica da Universidade do Porto and Spanish Ministerio de Ciencia e Innovacion under projects TIN2009-07519 and AIC10-D-000600.Vasconcelos, PB.; D'almeida, FD.; Román Moltó, JE. (2013). A Jacobi-Davidson type method with a correction equation tailored for integral operators. Numerical Algorithms. 64(1):85-103. doi:10.1007/s11075-012-9656-9S85103641Absil, P.A., Mahony, R., Sepulchre, R., Dooren, P.V.: A Grassmann–Rayleigh quotient iteration for computing invariant subspaces. SIAM Rev. 44(1), 57–73 (2002)Ahues, M., Largillier, A., Limaye, B.V.: Spectral Computations with Bounded Operators. Chapman and Hall, Boca Raton (2001)Ahues, M., d’Almeida, F.D., Largillier, A., Titaud, O., Vasconcelos, P.: An L 1 refined projection approximate solution of the radiation transfer equation in stellar atmospheres. J. Comput. Appl. Math. 140(1–2), 13–26 (2002)Ahues, M., d’Almeida, F.D., Largillier, A., Vasconcelos, P.B.: Defect correction for spectral computations for a singular integral operator. Commun. Pure Appl. Anal. 5(2), 241–250 (2006)Bai, Z., Demmel, J., Dongarra, J., Ruhe, A., van der Vorst, H. (eds.): Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide. Society for Industrial and Applied Mathematics, Philadelphia (2000)Balay, S., Buschelman, K., Eijkhout, V., Gropp, W.D., Kaushik, D., Knepley, M., McInnes, L.C., Smith, B.F., Zhang, H.: PETSc Users Manual. Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010)Chatelin, F.: Spectral Approximation of Linear Operators. SIAM, Philadelphia (2011)d’Almeida, F.D., Vasconcelos, P.B.: Convergence of multipower defect correction for spectral computations of integral operators. Appl. Math. Comput. 219(4), 1601–1606 (2012)Falgout, R.D., Yang, U.M.: Hypre: A library of high performance preconditioners. In: Sloot, P.M.A., Tan, C.J.K., Dongarra, J., Hoekstra, A.G. (eds.) Computational Science - ICCS 2002, International Conference, Amsterdam, The Netherlands, April 21–24, 2002. Proceedings, Part III, Lecture Notes in Computer Science, vol. 2331, pp. 632–641. Springer (2002)Henson, V.E., Yang, U.M.: BoomerAMG: A parallel algebraic multigrid solver and preconditioner. Appl. Numer. Math. 41(1), 155–177 (2002)Hernandez, V., Roman, J.E., Vidal, V.: SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems. ACM Trans. Math. Softw. 31(3), 351–362 (2005)Hernandez, V., Roman, J.E., Tomas, A., Vidal, V.: SLEPc Users Manual. Tech. Rep. DSIC-II/24/02 - Revision 3.1, D. Sistemas Informáticos y Computación, Universidad Politécnica de Valencia (2010)Saad, Y.: Iterative methods for sparse linear systems, 2nd edn. Society for Industrial and Applied Mathematics, Philadelphia (2003)Simoncini, V., Eldén, L.: Inexact Rayleigh quotient-type methods for eigenvalue computations. BIT 42(1), 159–182 (2002)Sleijpen, G.L.G., van der Vorst, H.A.: A Jacobi–Davidson iteration method for linear eigenvalue problems. SIAM Rev. 42(2), 267–293 (2000)Sorensen, D.C.: Implicit application of polynomial filters in a k-step Arnoldi method. SIAM J. Matrix Anal. Appl. 13, 357–385 (1992)Stewart, G.W.: A Krylov–Schur algorithm for large eigenproblems. SIAM J. Matrix Anal. Appl. 23(3), 601–614 (2001

Numerical Simulation of a Passive Control of the Flow Around an Aerofoil Using a Flexible, Self Adaptive Flaplet

© 2018 The Author(s) Self-activated feathers are used by almost all birds to adapt their wing characteristics to delay stall or to moderate its adverse effects (e.g., during landing or sudden increase in angle of attack due to gusts). Some of the feathers are believed to pop up as a consequence of flow separation and to interact with the flow and produce beneficial modifications of the unsteady vorticity field. The use of self adaptive flaplets in aircrafts, inspired by birds feathers, requires the understanding of the physical mechanisms leading to the mentioned aerodynamic benefits and the determination of the characteristics of optimal flaps including their size, positioning and ideal fabrication material. In this framework, this numerical study is divided in two parts. Firstly, in a simplified scenario, we determine the main characteristics that render a flap mounted on an aerofoil at high angle of attack able to deliver increased lift and improved aerodynamic efficiency, by varying its length, position and its natural frequency. Later on, a detailed direct numerical simulation analysis is used to understand the origin of the aerodynamic benefits introduced by the flaplet movement induced by the interaction with the flow field. The parametric study that has been carried out, reveals that an optimal flap can deliver a mean lift increase of about 20% on a NACA0020 aerofoil at an incidence of 20 o degrees. The results obtained from the direct numerical simulation of the flow field around the aerofoil equipped with the optimal flap at a chord Reynolds number of 2 × 10 4 shows that the flaplet movement is mainly induced by a cyclic passage of a large recirculation bubble on the aerofoil suction side. In turns, when the flap is pushed downward, the induced plane jet displaces the trailing edge vortices further downstream, away from the wing, moderating the downforce generated by those vortices and regularising the shedding cycle that appears to be much more organised when the optimal flaplet configuration is selected

Apoptotic cell-derived ICAM-3 promotes both macrophage chemoattraction to and tethering of apoptotic cells

A wide range of molecules acting as apoptotic cell-associated ligands, phagocyte-associated receptors or soluble bridging molecules have been implicated within the complex sequential processes that result in phagocytosis and degradation of apoptotic cells. Intercellular adhesion molecule 3 (ICAM-3, also known as CD50), a human leukocyte-restricted immunoglobulin super-family (IgSF) member, has previously been implicated in apoptotic cell clearance, although its precise role in the clearance process is ill defined. The main objective of this work is to further characterise the function of ICAM-3 in the removal of apoptotic cells. Using a range of novel anti-ICAM-3 monoclonal antibodies (mAbs), including one (MA4) that blocks apoptotic cell clearance by macrophages, alongside apoptotic human leukocytes that are normal or deficient for ICAM-3, we demonstrate that ICAM-3 promotes a domain 1–2-dependent tethering interaction with phagocytes. Furthermore, we demonstrate an apoptosis-associated reduction in ICAM-3 that results from release of ICAM-3 within microparticles that potently attract macrophages to apoptotic cells. Taken together, these data suggest that apoptotic cell-derived microparticles bearing ICAM-3 promote macrophage chemoattraction to sites of leukocyte cell death and that ICAM-3 mediates subsequent cell corpse tethering to macrophages. The defined function of ICAM-3 in these processes and profound defect in chemotaxis noted to ICAM-3-deficient microparticles suggest that ICAM-3 may be an important adhesion molecule involved in chemotaxis to apoptotic human leukocytes

Phytoplankton responses to marine climate change – an introduction

Phytoplankton are one of the key players in the ocean and contribute approximately 50% to global primary production. They serve as the basis for marine food webs, drive chemical composition of the global atmosphere and thereby climate. Seasonal environmental changes and nutrient availability naturally influence phytoplankton species composition. Since the industrial era, anthropogenic climatic influences have increased noticeably – also within the ocean. Our changing climate, however, affects the composition of phytoplankton species composition on a long-term basis and requires the organisms to adapt to this changing environment, influencing micronutrient bioavailability and other biogeochemical parameters. At the same time, phytoplankton themselves can influence the climate with their responses to environmental changes. Due to its key role, phytoplankton has been of interest in marine sciences for quite some time and there are several methodical approaches implemented in oceanographic sciences. There are ongoing attempts to improve predictions and to close gaps in the understanding of this sensitive ecological system and its responses